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I. INTRODUCTION 



In recent years, a series of remarkable results have been proven concerning the canonical quan- 
tization of linear scalar fields propagating on compact spatial manifolds, with a dynamics resem- 
bling that of a free field, but with an effective time-dependent mass. More precisely, the results 
apply to a field equation of the type 

X-Ax + s(t) X = 0, (1) 

defined on a static spacetime of the form Ex/, where / is an interval of the real line, £ is a compact 
Riemannian manifold of dimension d < 3, A is the Laplace-Beltrami (LB) operator on Z, and s(i) 
is a general function of time, subject only to rather mild conditions. For instance, for all purposes, 
it suffices that this function is twice differentiable, with a second derivative that is integrable in 
each compact subinterval of /. 

The results of Refs. yj-l3|] show that one can always find a Fock quantization for the system 
such that: i) the state that defines the Fock representation 1 is invariant under the isometries 2 of the 
spatial manifold S, and ii) the dynamics dictated by the field equation (QQ) is unitarily implemented. 

We notice that, although linear, the classical dynamics is nontrivial, owing to the presence of 
the time-dependent mass s(t). By unitary implementation of the dynamics we mean the following: 
given two arbitrary values of time, t and t', the linear symplectic transformation that corresponds 
to classical canonical evolution from t to t' is implemented in the Fock representation as a unitary 
operator. This is assumed to happen for all instants of time in (every connected component of) 
the interval /, with no further conditions. In particular, no continuity conditions are imposed, and 
therefore nothing is said about the existence of a well-defined Hamiltonian operator. 

In addition, and most importantly, the above Fock representation has been shown to be unique 
modulo unitary equivalence, in the sense that any other Fock representation with the same proper- 
ties of invariance under spatial symmetries and unitary dynamics is guaranteed to define a unitarily 
equivalent representation fllHJ. Note that these uniqueness results are obtained in the absence of 
time-translation invariance, which is a key ingredient in the standard uniqueness theorems, regard- 



1 This state is usually called the vacuum of the representation, although it does not necessarily correspond to an 
eigenstate of a Hamiltonian operator. 

2 More generally, one can consider the group of transformations that commute with the LB operator and are unitary 
in the Hilbert space of square integrable functions of the field configuration, with respect to the measure defined by 
the metric volume element. 
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ing the quantization of free fields in stationary (or static) spacetimes (see Refs. [|4j-l6|]). In me 
considered nonstationary settings, the mentioned new results single out a unique equivalence class 
of representations, ensuring the nonambiguity of the physical predictions, rather than selecting a 
definite Fock representation based on a specific vacuum. 

Though rather simple at first glance, these results find application in a variety of situations, 
including the quantization of inhomogeneous spacetimes, such as the Gowdy models (see Refs. 
|7-[l2]), the quantization of cosmological perturbations [Q, [13I 14], and the discussion of string 
dynamics in arbitrary plane wave backgrounds [15J. Another particular instance in which these 



results can be applied is the case of free fields in a nonstationary spacetime which is nevertheless 
conformal to a static universe, by means of a time dependent conformal factor. This is the situation 
found in Friedmann-Robertson- Walker (FRW) universes (with compact spatial sections), as well 
as in the very interesting case of the de Sitter spacetime 3 . In fact, in all of these examples, the 
use of conformal time combined with a suitable scaling of the original field variable (p transforms 
the original free field equation into and equation of the type CO)- The scaled field that satisfies this 
latter equation is x - Q l/2( f>, where Q is the (exclusively) time-dependent conformal factor. The 
new field possesses an effective time-dependent mass s(t) that depends on the conformal factor 
and on the mass of the original field. 

However, one can find claims in the recent literature stating that, for the case of the massless 
field in de Sitter spacetime, no scaling of the field variable allows for a unitary dynamics. That 
is the conclusion of Ref. QJ/ZU- One of the aims of the present work is to clarify this situation, 
showing with due care that those claims are unsound 4 . The standard conformal scaling of the field 
does indeed lead to a field formulation with the desired properties, allowing for a representation 
where one can reach a unitary implementation of the dynamics of the scaled field x, regardless 
of the value of the mass of the original field in de Sitter spacetime. In fact, in what concerns 
the possibility of a unitary implementation of the dynamics, the value of the mass parameter, 
positive, null, or even negative, is not relevant. The value of the mass is of course important for the 
existence of a fully de Sitter invariant Fock state of the Hadamard type. Actually, it is well known 
that the so called Bunch-Davies vacuum (or Euclidean vacuum), which is an 0(1, 4)-invariant 



3 For other criteria concerning free fields in 1 +1 dimensional de Sitter spacetime, see Ref. Ill 611 . 

4 The reasons can be traced back to an unsuitable choice of momentum field, as well as to the use of arguments based 
on the limit of infinite times, which is radically different from the limit in which the number of modes grows to 
infinity (for further details, see Sec. Pil l, 
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state, breaks down in the massless case H18H . Nonetheless, this does not prevent a consistent 
Fock quantization from being obtained, e.g. by means of 0(4), rather than 0(1,4), invariant 
states. We show explicitly that a Fock quantization can be achieved such that the dynamics of the 
massless field is unitarily implemented at the quantum level. This result is in full agreement with 
more general mathematical-physics results derived in the context of fields with a time-dependent 
mass. Taking advantage of such more general studies, one can moreover show that the obtained 
quantization is in fact unique, in the sense that any other 0(4)-invariant Fock representation which 
also allows for a unitary dynamics is necessarily unitarily equivalent. In particular, the Fock 
representation that naturally emerges from our approach is seen to be unitarily equivalent to the 
0(4)- invariant quantizations proposed some time ago by Allen and Folacci II 1911 . Moreover, by 
applying our approach to the free massive field case, we get a Fock representation which (i) admits 
a unitary implementation of the dynamics and (ii) has a vacuum state which is unitarily equivalent 
to the celebrated Bunch-Davies state. 

The paper is organized as follows. In Sec. HI] we introduce the classical setting for free, real 
scalar fields propagating in de Sitter spacetime. Conformal time is chosen as the evolution param- 
eter and a scaling of the original field variable by the spacetime conformal factor is performed to 
define the basic field variable. The issue of unitary dynamics for the massless field is discussed in 
Sec. Unl As the main result of the present work, a Fock representation permitting a unitary imple- 
mentation of the dynamics of the (scaled) massless scalar field is presented. Section |IV]is devoted 
to the study of the relationship between the so obtained Fock representation and those based on 
Hadamard states, namely the Bunch-Davies vacuum and the Allen-Folacci states. The relationship 
with another important set of states, namely the adiabatic states, is briefly addressed in Appendix 
lAl We summarize and discuss our conclusions in Sec. refsection: conclusions. Finally, Appendix 
|B] contains details of the calculations needed in the proofs of the unitary evolution and the unitary 
equivalence between different vacua. 

H. THE MODEL: FREE FIELDS IN DE SITTER SPACETIME 

The de Sitter spacetime is the maximally symmetric spacetime of positive constant curvature. 
It has the topology of R x S 3 and can be seen as the hyperboloid 

- x\ + x\ + x\ + x\ + x\ = H~ 2 (2) 
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embedded in flat five-dimensional spacetime (see e.g. Refs. Ill9l 2oL . The curvature of the space 



is r = 12H 2 . 

A system of coordinates (t, cr, 9, (p) in the whole space can be defined as follows: 

x = H~ l sinh(Ht), -oo < t < oo, (3) 

xy = H~ l cosh(Ht) cos(cr), < cr < n, (4) 

x 2 = H' 1 cosh(#?) sin(cr) cos(0), < 6 < n, (5) 

X?, = H~ l cosh(Ht) sin(cr) s'm(8) cos(^), < (f < 2n, (6) 

X4 = H~ l cosh(Ht) sin(cr) s'm(8) sin(^). (7) 
In these coordinates the metric reads 

ds 1 = -dt 1 + H~ 2 cosh 2 (Ht){dcr 2 + sm 2 (cr)[de 2 + sin 2 (6>) dy 2 ]}. (8) 

The de Sitter spacetime is conformal to the static universe R x 5 3 . To see this, let us introduce the 
conformal time 

T) = 2 arctan(e m ), r\ e (0, n). (9) 
In the new coordinate system, the metric takes the form 

ds 2 = a 2 (Tj){-drf + da 2 + sin 2 (cr)[^ 2 + sin 2 (0)V]}, (10) 

where one recognizes the metric of the static universe and the time-dependent conformal factor 

a{r 1 ) = -±—. (11) 
H sin T] 

Let us consider the propagation of a free (minimally coupled) real scalar field (f> with mass m 
and dynamical equation 

(□ - m 2 )(f> = 0, (12) 

where □ stands for the d'Alambertian associated to the spacetime metric. In the coordinates 
(t, cr, 6, if), this field equation becomes 

d 2 (p + 3H tmh(Ht)d t (f) + l~ + m 2 ^ = 0. (13) 

Here, A is the LB operator on 5 3 . We now change to conformal time and introduce the scaled 
field 5 

X = a<p. (14) 



5 Similar scalings can also be performed in less than three spatial dimensions to reach a field which admits a unitary 
dynamics fl. 



We then obtain the new field equation 

X+[-A + (m 2 -2H 2 )a 2 + l]x = 0. 



(15) 



The dot stands for the derivative with respect to the conformal time 77. Note the absence of terms 
containing the first time derivative of the field. 

The Lagrangian density corresponding to the field x ( U P t0 tota l ti me derivatives) is 



L-i 

2 



The canonical momentum conjugate to^ can then be taken as 



(16) 



P X =X- 



(17) 



Alternatively, one might have started from the canonical pair (0, = a 3 d t (p) and performed the 
time-dependent canonical transformation 

P, 



X = acf>, P x h a(p. 



(18) 



Let us return to the field equation (fT3T > and decompose the field x m eigenmodes of the LB 
operator: 



k,t,m 



where Y ktm are the 5 3 -spherical harmonics, which satisfy the eigenvalue equation 



&Y ktm = -k(k + 2)Y Um , 



(20) 



and provide an orthonormal basis for the space of square integrable functions on S 3 . In these 
formulas, the integer k takes values from to 00, i varies from to k, and m ranges from —I to I 
(see e.g. Refs. i2lLl22|l ). In the following, we will use the notation q k to collectively denote all the 
modes qu m corresponding to the same value of k. 

Introducing the decomposition (fT9l) in the field equation (fT5l) and using the orthogonality rela- 
tions for the harmonics, we obtain the dynamical equation of each mode: 

m' 2 -2 



(k + \y + 



sin 2 1] 



<?k = 0, 



(21) 



where m' = m/H. These equations are the same for all the values of t and m that share the same 
value of k. Note also that, from P x = x, it follows that the momentum canonically conjugate to 
the mode variable q k satisfies the equation p k = q k . 



The general solution to the equations of motion (|2TT) is well known. In fact, a change of variable 
^k(^) = (sin rf) ll2 f{- cos 77) transforms equation (|2TT > into 



(>-.-) 



£L_ 2 iL + 

dy 2 dy 



— — m 



f(y) = o, 



(22) 

1/2 



4 l-y 2 \4 

with y = - cos 77. This is a Legendre equation of degree v = k + 1/2 and order // = (9/4 - m' 2 ) ' , 
whose independent solutions are the associated Legendre functions Py and <2^ (see e.g. Refs 



24j]). The general solution to the equations of motion (|2T|) is therefore 
q\M) = A k Vsin^(- cos 77) + 5 k ^%mrjQ^{- cos 77), 



(23) 



where Ak and 5t are arbitrary (complex) constants. 

In the canonical formalism, the general solution to the corresponding Hamiltonian equations 

4k = Pk, Pk = ~ [(k + l) 2 + (m' 2 - 2) sin" 2 77] q k , 

can then be written in the form 

Ry(- cost/) 5y(-cos?7) 
^(-cos?7) S^- cos 77) j 



(24) 



<7k 




V 




= M k (T]) 




M,(/7) = 


,/>k v 







(25) 



where 7?y and Sy are given by 



R^(-cost])= -y'sin r]P^(- cos 77), S^(- cos 77) = V sm7 72v( - cos 7 7)- 
By using the relation 

^ = ^ T «(y)-(v + ^_ 1 (y)L 



(26) 



(27) 



valid also for Q^(y) (see e.g. Refs. |23j, |24J]), i?y(- cos rj) and 5y(-cos77) can be expressed as 
follows 

^(-cos/7) = — =f(v+ l/2)cos77^(-cos77) + (v+^_ 1 (-cos7 7 )l, (28) 
y sin 77 

^(-cos/7) = ^^[(v+l/2)cos77^(-cos ? 7) + (y+//)^;_ 1 (-cos77)l. (29) 
y sin 77 

Note that one can write the above matrix elements directly in terms of the argument cos 77, 
instead of - cos 77, since [24] 



/*(-*) = cos[(v + iu)n]P^x)--sm[(v + iu)n]Q'i(x), 

n 



n 



Q>X-x) = -cos[(v + M)nW v (x)--sm[(v + iu)n]P^x). 



(30) 
(3D 
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It is not difficult to see that the same expressions hold for Ry and S f [replacing P\,(+x) and Q^{+x) 
with Ry(+x) and Sy(+x), respectively]. In particular, for the massless case (// = 3/2), we get 
simply 

R^(-x) = (-l)^(x), Sft-x) = (-l)* +1 S£(x), (32) 

where we have used that v = k + 1/2. 

Time evolution from an arbitrary reference time 770 to any another time 77 is then given by the 
canonical transformation 



, Pk(r]) , 



= [4(770, 77) 



/ <?k(?7o) 
£>k(?7o) , 



(33) 



where U k (rj , 77) = M k (r])M k l (r] Q ). 



HI. QUANTIZATION WITH UNITARY DYNAMICS 

We will now show that one can find a Fock quantization which allows a unitary implementa- 
tion of the dynamics of the field x, i.e., a unitary implementation of all transformations (1331 , for 
arbitrary values of 77 and 770. 

A Fock quantization is defined by a choice of a complex structure on phase space 6 , which 
is tantamount to a choice of sets of creation and annihilation variables (up to irrelevant changes 
which do not mix both sets). Let us introduce the classical (complex) variables 

a k = - — {co k q k + ip k ) , a k = * — (co k q k - ip k ) , (34) 
y2io k \2u> k 

where the frequency co k is chosen to match the time-independent part in the equations of motion 
(T2T)) . i.e. co k = k + 1. Naturally, these variables satisfy canonical Poisson brackets {ak, a* k ,} = i6 kk >. 

If we now declare that the variables a k and a k are to be quantized as creation and annihilation 
operators, then we single out a particular Fock quantization. In other words, the complex structure 
J that determines our particular Fock quantization is defined by J(a k ) = ia k , J(a* k ) = -ia* k (see 
Ref. [2] for details on the Fock quantization). 



6 Remember that a complex structure is a map on phase space that preserves the canonical structure and whose square 
is minus the identity. For the construction of a Fock representation, one demands that the complex and the canonical 
structures be compatible in the sense that a suitable composition of their actions be positive definite. 
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Classical time evolution, written in terms of the variables a^, at, is given by 



a k (r]) 


= WkOjo, J]) 








, afro) , 



T k U k (T] ,j])T k l 



(35) 



where 



f \ 
a>„ i 

\0) n -i 



(36) 



is the matrix corresponding to the change of variables ([34V Since the transformations (|35T ) are 
canonical, they necessarily take the form of a general Bogoliubov transformation, i.e., they can be 
written in the form 



' a k (f] , rj) p k (i] , rj) 
/3* k (T]Q,T]) a* k (rj ,r]) j 



with 



\a k (j]o,n)\ ~ \Mm,ri)\ = 1, 



(37) 



(38) 



independently of the particular values of k, rj , and rj. 



Standard results Q25 



2611 now show (see Ref. [|2J] for details) that the dynamics is unitarily 
implementable (in the above Fock quantization) if and only if the functions fi k in Eq. (1371) are 
square summable; that is, if and only if 



e 



TjTjTj ^(Vo,ri)\ 2 = J](k + l) 2 ^(77o,/7)l 2 < °°, 



(39) 



k=Q {=0 m=~C k=0 
\2 



where the degeneracy factor (k + 1) counts the number of degrees of freedom with the same 
dynamics. The fulfillment of this summability condition depends on the asymptotic behavior of 
the functions /3 k (rj , rj), for large values of k. This in turn depends on the asymptotics of Ry and 
Sy dictated by the ultraviolet (large values of the degree v = k + 1/2) behavior of the Legendre 
functions Py and <2v- The expansions of Py and (£, are as follows I124ll 

^T(v + + ^ rOu + y -h 1/2) cos f(v + 7 + 1/2)77+ |(2j - 1) + f 1 



P$(COSTj) = \ — 



/sin rj 



TQi-j+ 1/2) F(j + 1) r(v + j + 3/2) (2 sin rjV 



ee(cos77)= \ — 



2 r(v + fx + 1) y rQu + j + 1/2) cos [(v + j + 1/2) r, - | (2j - 1) + f ] 
Zj (_ 1 y TOu - j + 1 /2) fC/ + 1 ) Hv + j + 3/2) (2 sin rjV ' 



Vsin?7 ; . =0 

These formulas provide asymptotic expansions (up to arbitrary order) valid not only for real values 
of the parameter /i, but also for complex ones. Thus, the range of validity of the asymptotic 



expansions covers the whole set of possible values of m, including, of course, the massless case 
0/ = 3/2). 

Let us momentarily focus on the massless case. The asymptotic behavior of the Bogoliubov 
coefficients /^(t/o, rf) can now be deduced, taking into account the matrices M k , relations (12~8l [29}, 
the matrices £4 an d T k , identities 021 , and the asymptotics for the functions Pl /2 and Ql /2 . A 
lengthy but straightforward computation, detailed in Appendix [Bl shows that fik(rjo, rf) is of order 
O (k~ 2 ^), for large k, i.e., the asymptotic behavior when k — » oo is given by 

p k ( m ,rj) = 0{k- 2 ), Vt7,77o. (42) 

Thus, it follows that the summability condition (1391) is satisfied for all values of -qo and 77, and there- 
fore the dynamics is unitarily implemented in the considered Fock representation 7 . This result is 
in complete agreement with the general results proven in Refs. , and disproves the conclusion 
of Ref. jn\\ . where it is claimed that one cannot attain (by means of a Fock quantization) quantum 
unitarity of the evolution for the massless field in de Sitter spacetime, independently of the field 
redefinition <p —* fif)<f>. 



The calculations in Ref. [17] are based on a specific choice of momentum for the scaled field: 
the momentum obtained just by the inverse scaling. Nevertheless, one can also adopt other choices 
of momentum in order to obtain a canonical pair, while respecting the linearity of the system. 
Namely, one can allow for a time-dependent linear contribution of the field configuration to the 
momentum. This kind of time-dependent linear transformations affect the dynamics of the basic 
field variables. Actually, time-dependent canonical transformations of this type have been shown 
to be crucial to arrive at a unitary evolution [12711 . More specifically, a unique transformation is 
admissible when the dimension of the spatial hypersurfaces is greater than one, as it is the case 
here. The appropriate canonical transformation (including the field scaling) that leads to unitary 
evolution is that provided in Eq. (fTHT ). In view of the transformation, it is particularly convenient 
to describe the system in conformal time, because then the privileged momentum is given just by 
the time derivative of the scaled field and the field equations simplify considerably, reflecting the 
conformal nature of the spacetime metric. In Ref. [17], nonetheless, this conformal nature was 
not explored. Only in the concluding section of that work the conformal time was considered, 



7 Since the complex structure J depends on the LB operator only, which in turn is an <9(4)-invariant object, our unitary 
Fock representation is, in addition, an (9(4)-invariant quantum theory. Let us note also that the same asymptotic 
behavior and therefore the same conclusions apply for any other value of the mass (see Appendix IBl. 
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presenting just a heuristic argument to support that the conclusions about the nonunitarity of the 
dynamics were valid as well for conformal time. That argument is however not correct. The 
argument goes as follows. When the scaling ([141) and the conformal time r\ are used, the field 
equation that one obtains for the massless case, namely Eq. (fT51) with m = 0, is of the form of a 
field with a time-dependent mass which is always negative. That is indeed the case, since we get 
from Eq. (fBT) 



1 - 



sin 2 77 



X = 0. (43) 



The time-dependent mass m(rf) = 1 - 2/ sin 77 is not only strictly negative, but moreover blows 
up when 77 — > 0, 7r (which corresponds to t — » +00). It was then argued in Ref. Jn| that, given 
that the time-dependent squared frequencies of the harmonic modes, namely (k + l) 2 - 2/ sin 2 77, 
are all negative in the limit t — > +00, this would introduce a non-oscillatory behavior in that limit 
that would cause the failure of unitarity. However, the above argument does not really pose any 
obstruction to the unitary implementation of the dynamics. It rather points out that the ultraviolet 
limit, in which the infinite number of modes of the system come into the scene, and the limit of 
infinite time t are radically different. In fact, unitary dynamics means the unitary implementation 
of all the evolution transformations (1331 between any two finite values of time. Actually, when 
one considers the evolution between two instants of time, 770 and 77, one should look to the values 
of (k + l) 2 - 2/ sin 2 77 in the limit of large k, and not in the limit of large t, because the dynamical 
transformation is sensitive only to the values in the interval [770, 77], since the equation of motion is 
local in time. What happens is that for all k greater than some (maybe large, but finite) order ko 
(which depends on 770 and 77), the values of the squared frequencies (k+ 1) 2 -2/ sin 2 77 are positive in 
all the considered time interval. It is true that, for a finite number of modes, the evolution (between 
finite times) takes place with a negative time-dependent frequency. But this does not affect the 
possible unitary implementation of the dynamics in a fundamental way, because linear dynamics 
in finite dimensions is always unitary, as it is granted by the Stone-von Neumann uniqueness 
theorem J28]. 

For a full clarification of this situation, let us consider the case where the time dependent term 
2/ sin 2 77 in the field equation is replaced simply by 2/?7 2 . This corresponds of course to the limit 
t — » +00 in the scale factor (fTTT ) and is therefore actually physically relevant. Moreover, the field 
equation obtained with this replacement in the massless case, namely 



-7 



X = 0, (44) 
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keeps precisely the qualitative features of the original equation (1431) that are involved in the argu- 
ment sketched above about the behavior of the mode frequencies in the region of large \t\. The 
advantage of Eq. (1441) is that it can be solved explicitly in terms of elementary functions, thus 
making the whole discussion fully transparent. 

In fact, the equations of motion that we obtain for the harmonic modes, corresponding to Eq. 
(|2TI) (with m = 0) are now 



2 

n 

One can readily check that the general solution is of the form 



qy = 0. (45) 



qm = Ak ( !ifi±lM _ cos[( * + ,)„]) + Bt L *°f ±M _ sin[( * + m ) , (46 ) 

where A k and 5 k are arbitrary complex constants. One can now simply follow the procedure de- 
scribed above and obtain the corresponding coefficients /^(t/o, r/) in the Bogoliubov transformation 
(1371) . With A77 standing for 77 - 770, we get in this case 

AC,., n) - ^ (1 - 1 + j^^) com * DA,] 



The leading behavior is thus of order o(k~ 2 ^j, and it follows that the summability condition 
(l39l) is indeed satisfied. This was the conclusion expected from more general results concerning 
the unitary implementation of time evolution fllsl, and confirms, in particular, that the sign of the 
mass term is irrelevant 8 . 

We have therefore constructed an 0(4)-invariant Fock representation for the massless free scalar 
field in de Sitter spacetime, with unitary dynamics. In what follows, we will investigate the rela- 
tionship between our Fock vacuum and the family of 0(4)-invariant Hadamard states characterized 
by Allen and Folacci Before doing that, let us stress that, although we have focused on the 
massless case, the unitarity result holds as well for any massive free field; i.e., the Fock repre- 
sentation, defined by the choice (1341) of annihilation and creation operators, provides a quantum 
description where time evolution of the free massive field admits a unitary implementation. That 
this is so can be directly verified from the asymptotic behaviors of Py and which, for any (con- 
stant) value of the parameter jx (including complex numbers), provide a beta coefficient satisfying 



A comment in this respect already appeared in Ref. |15 1 
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Eq. (1421) . Hence, given a scalar field with m > (what is more, with any real -even negative- 
value of m), there exists at least one Fock representation where time evolution is implemented as 
a unitary operator. Further details of this proof can be found in Appendix El 



IV. EQUIVALENCE WITH THE ALLEN-FOLACCI'S <9(4)-INVARIANT STATE 



An alternate way of defining a Fock quantization consists in selecting a particular set of com- 
plex mode solutions {u^} to the equations of motion (12TT) . These solutions are normalized so that 
they satisfy the condition 

u k u k - u* k iik = i (48) 

on a given Cauchy surface, say rj = rj . The symbol * denotes complex conjugation. Once such a 
set of solutions is chosen, one can write the general solution to the field equation (fTBI) as 



(49) 



The Fock quantization is then performed by declaring that b k and b* k are to be quantized as the 
annihilation and creation operators of the Fock representation (see Ref. Ii29i1 for details and a nice 
account on Fock states on homogeneous and isotropic spaces). In this description, the Fock quan- 
tization presented in the previous section corresponds to the choice of mode solutions determined 
by the following initial data: 



1 UOk 

u k (t] ) = — — , u k (r] ) = -i J — . 
y2a> k V I 



(50) 



For massive free fields in de Sitter spacetime, the Fock quantization is usually carried out using 
the unique solution which is de Sitter invariant, i.e., invariant under the full 0(1,4) group, and 
satisfies the Hadamard criterion. The corresponding set of mode solutions is of the form (|23l) . with 

HQ 

Bk = A iAk , Ak knk- M+ 3/2) 

n ^4T(k+iu + 3/2) V J 

The vacuum of the corresponding Fock representation is known as the Bunch-Davies, or Eu- 
clidean vacuum. Explicitly, the mode solutions determining this Euclidean quantization can be 
written as 



R^{- cos 77) iS^{- cos 77) 

n 



A k e 



i(v+fi)n 



R^icosrj) + -iSt(C0ST]) 
n 



(52) 
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where in the last equality we have taken into account Eq. (1261) as well as relations (13013 11 1. Given 
a Cauchy surface, say rj = rj , the momentum canonically conjugate to ^(770) is 



2 . 

i^(cos 77) + -iS%cos rf) 
n 



In the case of zero mass, it is a well known fact that the Euclidean vacuum breaks down, i.e., 



(53) 

'7(1 



there is no de Sitter invariant Hadamard vacuum in this case [18]. As explained in Refs. [18, 
this is due to the dynamics of the zero mode only. One can in fact easily check that taking m = 0, 
and hence n = 3/2, in the above expressions (15 11521 ). one gets perfectly well-defined solutions 
for k ^ 0, whereas the corresponding expression for k = becomes meaningless. To obtain a 
complete set of well-defined solutions, and therefore a well-defined quantization, one only needs to 
derive proper solutions for the zero mode [or just quantize the zero mode in a different, consistent 
alternate way] . Independent solutions to the zero mode equation of motion with m = are 1 / sin rj 



and (77/ sin?/) - cost/. It is shown in Ref. 111911 that one can arrive at a one-parameter family of 
solutions for the zero mode such that [together with the solutions (|52l) for k ^ 0] 0(4)-invariant 
Hadamard vacua are obtained. 

The question naturally arises of whether or not the quantization presented in the previous sec- 
tion is unitarily equivalent to those corresponding to the 0(4)-invariant states constructed by Allen 
and Folacci. Note that the unitary equivalence between those representations depends only on the 
behavior of the states for large values of k, and therefore the particular quantization used for the 
zero mode is irrelevant (provided of course that we are using for the zero mode a standard quanti- 
zation that satisfies the Stone-von Neumann conditions, and not e.g. a polymer type quantization). 

Given two sets of mode solutions [^ k } and {x' k }, determined by initial conditions ICtk^klL an d 
{Otk'A'UL' the two corresponding Fock representations are unitarily equivalent if and only if the 



following set of /3k coefficients is square summable (see e.g. Ref. [29]): 



fik = i |> k (?7okk(?7o) - Xk(m)X k (m)] ■ (54) 

Again, for the sets of solutions that we are considering, the coefficients /3k depend only on the 
index k, and not on the full set of labels k. So, the summability condition is still of the type (|39l , 
with a degeneracy factor (k + l) 2 . 

Let us now fix a Cauchy surface, r\ = 770, and evaluate the coefficients /3 k relating our 
data/representation (l50l ) with the Allen-Folacci data/representation (Xk,Xk)\ri Q with ju = 3/2 [see 
Eqs. (1521531) 1. Using relations (1281291) and the asymptotics for the functions Pl /2 and Ql /2 , one can 
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check (see Appendix |B] for further details) that k in Eq. (|54l) is of order O (k~ 2 ) in the ultraviolet 
regime k — » oo (or, equivalently, when v — > oo). This asymptotic behavior is sufficient to satisfy 
the summability condition, and it is therefore proven that the Fock representation discussed in the 
previous section is unitarily equivalent to the representation defined by the Allen-Folacci vacua. 

In addition to the massless case, which corresponds to yu = 3/2, one can further check (as we do 
in Appendix |B]) that the asymptotic behavior /? k ~ O (^~ 2 ) still holds for any m > 0, case in which 
the solution (|52|) corresponds to the celebrated Bunch-Davies (or Euclidean) vacuum. Hence, the 
Fock representation of Sec. |III] is unitarily equivalent to the representation based on the Bunch- 
Davies vacuum. Let us again remark that the quantization of Sec. Hill supports a unitary dynamics 
for massive fields as well. We have thus proven, also for massive fields, that the quantization based 
on the Bunch-Davies vacuum, which follows from the requirement of full de Sitter invariance and 
the Hadamard condition, is unitarily equivalent to the quantization obtained from the requirement 
of unitary implementation of the dynamics. 

The fact that both viewpoints -the one using the Hadamard condition and the one based on 
unitary dynamics- select the same equivalence class of representations [when the scaled fields is 
used] can only be considered as a reassuring result, connecting two a priori distinct approaches. In 
the unitary dynamics perspective, one imposes only the existence of unitary transformations that 
implement the classical time evolution between any two (regular) instants separated by a finite (as 
opposed to infinitesimal) interval of time, with no extra requirements, like e.g. continuity or any 
pre-assigned local form of the vacuum. When one uses the Hadamard condition, an apparently 
stronger condition is imposed, that fixes the local singularity structure of the vacuum state. As it 
is well known, this condition is strong enough to guarantee the regularization of the stress-energy 
tensor, which was in fact the original motivation to adhere to it. It does not seem at all clear 
that those two viewpoints should lead to fully equivalent quantizations, and the fact that they do 
constitutes an interesting result by itself. 

V. CONCLUSIONS AND DISCUSSION 

We have explicitly shown that there exists a Fock quantization of the massless scalar field in de 
Sitter spacetime admitting a unitary implementation of the time evolution. Like in other situations 
considered in the literature, this involves a scaling of the field variable, using the conformal factor, 
and the introduction of a suitable momentum field, given by the conformal time derivative of the 
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scaled field. Our result disproves previous statements in the literature (11711 claiming precisely that 
it is impossible to attain a Fock quantization with unitary dynamics in this case, by means of the 
scaling that has proven to be so effective in other situations fllml uMllll. In addition to providing a 
direct proof of the unitary implementation of the dynamics of the massless scalar field in de Sitter 
spacetime, we have analyzed a completely solvable toy model, which further clarifies the viability 
of a unitary evolution. 

Besides, we have seen that the same Fock representation also supports a unitary implementation 
of the dynamics of massive free fields. It is worth remarking that general results allow us to ensure 
that the Fock representation depicted in Sec. |III] is the unique one (modulo unitary equivalence) 
admitting a unitary implementation of the time evolution of free fields in de Sitter spacetime. For- 
tunately, there is no tension between this result, on the one hand, and the uniqueness provided 
by imposing the Hadamard criterion, on the other hand. Let us recall that, for free scalar fields 
propagating on spatially compact universes, the Hadamard approach selects a unique preferred 
representation of the canonical commutation relations D5|]. Thus, in particular, a Klein-Gordon 
field (p in de Sitter spacetime has a unique quantum Hadamard representation. More specifically, 
given a massive free field in de Sitter spacetime, there exists a unique 0(1, 4) -invariant Fock vac- 
uum state satisfying the Hadamard condition: the Bunch-Davies (or Euclidean) vacuum state. In 
the massless case, there are instead infinitely many 0(4)-invariant Hadamard vacua, differing in 
their particular zero mode sector parametrization [11911 . However, since the discrepancy between 
vacua involve just a finite number of degrees of freedom, the Stone-von Neumann uniqueness 
theorem guarantees that the family of 0(4)-invariant Hadamard vacua is, in fact, a family of uni- 
tarily equivalent states. Thus, in de Sitter spacetime, there is a unique 0(1, 4)-invariant Hadamard 
quantization of a massive Klein-Gordon field, and a unique (equivalence class of) 0(4)-invariant 
Hadamard quantization(s) of the massless Klein-Gordon field. Under the time dependent canon- 
ical transformation (fTST ). the unique Hadamard quantization of the (m > 0) Klein-Gordon field 
cp determines, in the scaled field description a quantum theory which is characterized by the 
Cauchy data (|52I53I) . Our results of Sec. |lV]show that for any (nonnegative) value of the mass 
parameter, this translated (unique) Hadamard quantum theory defines a Fock quantization which 
is unitarily equivalent to our Fock representation; i.e., the so translated Hadamard quantization and 
the Fock representation with unitary dynamics provide exactly the same physical predictions. This 
equivalence eliminates any concern about a possible tension between the requirement of a unitary 
time evolution (together with the invariance under the spatial symmetries) and the Hadamard con- 
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dition in order to select a Fock representation. 
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Appendix A: Equivalence with adiabatic states 

For the sake of completeness, we devote this appendix to discuss the relationship between our 
Fock quantization with unitary dynamics and the use of adiabatic states in de Sitter spacetime. As 



it is known, adiabatic states were introduced in the late sixties by Parker [30] to bring forward the 
best possible definition of particles in expanding universes scenarios. In the framework of closed 
universes (specifically, in FRW cosmologies with k = +1) two important mathematical-physics 
results exist concerning adiabatic states: (i) the family of adiabatic vacua is a set of unitarily 



equivalent states |29J], and (ii) adiabatic states are unitarily equivalent to Hadamard states (for 
m > 0) [31]. We thus have, in particular, that the (Hadamard) Bunch-Davies vacuum state is 
unitarily equivalent to an adiabatic state. Given the equivalence between the Bunch-Davies state 
and the unitary Fock vacuum state (namely, the vacuum of our Fock quantization with unitary time 
evolution), established in Sec. [IV] we then conclude that our Fock vacuum is unitarily equivalent 
to an adiabatic state. We will explicitly show here that our Fock state is unitarily equivalent to the 
zeroth order adiabatic vacuum state. 

In order to properly compare the zeroth order adiabatic vacuum state, defined in the (f>- 
description, with the unitary Fock vacuum state, defined in the scaled ^-description, we will 
proceed in three steps. First, we will choose x as me basic field variable. Next, we will trans- 
late the Cauchy data defining the zeroth order adiabatic state to the ^-description, and finally we 
will compare the result with our unitary Fock vacuum state, defined by the Cauchy data (|50l) . 

Let us start by considering the field equation (fT3l for a free scalar field propagating in de Sitter 
spacetime. By performing a mode decomposition of the field, we get that the time-dependent part 
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of the mode solutions, Vk, obey the second-order differential equation 

dW + 3 (^) d t v k + w 2 k v k = 0, wl = k -^± + m\ (Al) 

The modes Vk must satisfy, in addition, the normalization condition 

Qy.Pl ~ QiP* = i, Vk, (A2) 

where Q k = v k (t ) and P k = a 3 d t v k (t ) are, respectively, the configuration and momentum coeffi- 
cients of the field </> on the Cauchy surface t = to. 

In order to introduce adiabatic vacuum states, let us consider solutions to Eq. (IA1I) of the form 

v k (t) = —^=exp(-i f e k (t)dt\ (A3) 
where k are real positive functions which, according to Eq. (IA11 ) and Eq. (IA3I) . must satisfy 



, 3(d t a\ 2 2>d 2 a 3 /d f k \ 2 \d 2 t ® k 

This equation can be solved via an iterative process whenever a finite time interval and a suf- 
ficiently large k are considered [29]. Thus, starting the process with k O) = w^, we get in the 
left-hand side of Eq. (|A4|) the (r + l)-th function k r+1) by plugging in the right-hand side the 
preceding r-th solution Q k \ 

An adiabatic vacuum state of r-th order is a Fock state constructed from a solution {vk} to Eq. 
(lATT) with initial conditions (v k \t Q ), d t v k r \t )), where 9 

v k r \t ) = 1 exp Li f ° e ( k r \t)dt\ (A5) 

We have k O) = w k = [k(k + 2) + (ma) 2 ] 1 /a for the zeroth order, so that the Cauchy data of 
the corresponding adiabatic vacuum state is given by 

„2 



(Jk - v k , F k - -a v k 



YYl 

1 H ) d t a + iawu 

2w\ 



(A6) 



9 It is worth remarking that adiabatic vacuum states are independent of the values chosen for t and to. Indeed, different 
choices of t in Eq. ( IA5t just introduce irrelevant phases, whereas distinct elections of the reference (initial) time fo 
lead to equivalent vacuum states J29I . 
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Now, we translate the Cauchy data (IA6I ) to the ^-description via the time-dependent canonical 
transformation (fTST ). 



(0) 

q k = aV k 



1 



(0) 



'70 



= -xv k 



m , 
a— + 2za w# 



'70 



(A7) 



where t]q = 2 arctan(e ff '°) [see Eq. ©]. 

Recall that the Cauchy data defining the unitary Fock vacuum state are given by Eq. (|50t . 
(uk(no) = (2w /t )~ 1/2 , UkiVo) = -'(^/ £ /2) 1/2 ). Thus, the antilinear part of the Bogoliubov transforma- 
tion relating the Cauchy data (|50l ) and (IA7I ) reads 



1 



((o k q k - ip k ) 



2(1 -x 2 ) 1 / 4 



. d t a 3 



(A8) 



V2o^ -v- 
where Xk = m/Wk- The unitary Fock vacuum and the adiabatic vacuum will be unitarily equivalent 
states if and only if fi k is square summable. It is a simple exercise to see that, in the asymptotic 
regime, which corresponds to Xk « 1, the ultraviolet behavior of beta is /?k ~ 0(k~ 2 ). As a 
consequence, we conclude that the unitary Fock vacuum state is equivalent to the zeroth order 
adiabatic vacuum. Now, since in closed FRW spacetimes any two adiabatic states of distinct order 
are unitarily equivalent, we have in fact that our equivalence result extends to adiabatic states of 
arbitrary order. 



Appendix B: Ultraviolet behavior of the beta functions and coefficients 

In this appendix, we will detail the derivation of the ultraviolet behavior, i.e., the behavior for 
large k, of the beta functions fi k (rio, rf) of the evolution transformation, defined in Sec. Hill and 
of the beta coefficients 0^ of the canonical transformation that relates the vacuum selected by the 
unitary evolution criterion and the Allen-Folacci vacuum, discussed in Sec. [IV] 

Let us start by detailing the behavior for large k of the beta functions, i.e., the antilinear part 
of the evolution transformation f Uk(j] ,rj) introduced in Eq. (|33T > . We first recall the expression 
of this transformation on the creation and annihilation variables in terms of the evolution of the 
configuration and momentum modes, 

17) = T k Uk(7]o, n)^ 1 = TkM k (rj)M- k \ nQ )T^ . (Bl) 

The matrices M k and T k are given in Eq. (1251) and Eq. (l36l) . respectively. It follows from this 
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expression that the beta functions are 

r(^)r(^±i) 



2^ +1 r(^)r(^) 

where we have introduce the notation 



A DM CM 

AE$$t - AS ^ + i— ^ + ia> k ASW 



(B2) 



AFtG* = F?(- cos nY%(r cos 770) - G£(- cos iflFft- cos 770), (B3) 

for any function Fy and Gy in {F£,Sy,Fy,Sy}. Besides, to compute the inverse matrix M k l , we 
have employed the determinant of M k : 

det M k (rj) = sin 2 n W{P*(- cos 77), Q?(- cos 77)} = (B4) 

Here, f W{-, •} denotes the Wronskian. By considering the relations between the functions Py(x) 
and 2t'(x) with F^(-x) and <2^(-jc), given in expressions (1301311 ). one gets 

AR*$$ = S?(cos 77)^(cos 770) - fl£(cos 77)5^(008 770), (B5) 

and similarly for A5%, ARySy, and AS%. 

Using then Eqs. (14014 II) for the asymptotic expansion of the associated Legendre functions, as 
well as the definition (l26l) of the functions Ry and Sy, and the expressions (1281291 ) of their time 
derivatives, one finds that the beta functions have the asymptotic behavior 

*{ ^ r (~ 4~ ) r H i~ ) r 2 (fe + J u + 3/2) Vow 

*>> = - 2 ^ lr (2^ )r (2^r(fc + i)r(ifc + 2) £ ^ UoX (B6) 

where the functions 33 ; (?7, 770) are of asymptotic order O y^y Employing the Stirling formula for 



the asymptotic behavior of the Gamma function [23], 



T(z+ 1) ~ ^Jlnzizlef, (B7) 



it follows that the time-independent coefficient that multiplies the functions £8/(77, 770) in Eq. (1B6I) 
behaves asymptotically as 0(1). Therefore, the beta functions have the same ultraviolet behavior as 
the first nonvanishing function £8^(77, rj ). We will explicitly show that both 23 (t7, 770) and £81(77, 770) 
vanish, proving that the beta function is (at least) of the asymptotic order O (^ 2 ) ~ O (k~ 2 \ 

The function 93 (t7, 770) can be deduced from Eq. (1B2I) by considering only the leading asymp- 
totic contributions of the terms of the form AFyGy in that expression. For this, one can use the 
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asymptotic expansions of Ry, S f, and their time derivatives, which can be obtained from Eqs. (1401) 
and (|4TT) as well as Eqs. (1281) and (|29l) . With a bit of calculus, and using trigonometric relations, 
one gets the following four leading contributions, which annihilate each other in pairs: 

m) = ~ sin i(k + l)(*7+»7o) + /«r] + sin [(k + l)(i]+T]o) + fur] 

+/cos [[(k + l)(j]+Tj ) + fin] - / cos [[(k + 1)(ji+tj ) + fin] = 0. (B8) 

On the other hand, the expression of the function 2?i(?7, 770) can be computed in the same way, 
but considering now the next leading contributions of the terms of the type AFyGy. These are 
obtained from the products of the leading and first-subleading contributions in the asymptotic ex- 
pansions of Ry,Sy, and of their time derivatives. As before, making use of trigonometric relations, 
one obtains 

8+ _2 % i( - ) = _L { cos [{k + 2)7/ + (k + 1)770 + fur] - cos [(k + 2)77 + (k+ 1)770 + jot] 
4/r - 1 sin 77 

- i sin [(/c + 2)77 + (k+ 1)770 + jU7r] + i sin [(/c + 2)77 + (k + 1)770 + //7r]} 

+ — - — {cos Uk + 1)77 + (k + 2)770 + fin] - cos [(/c + 1)77 + (k + 2)r/ + fin] 
sin 770 

- i sin [(k + 1)77 + (k + 2)770 + fin] + i sin [(/c + 1)77 + (k + 2)rj + fin]} = 0. (B9) 

Therefore, both 23 (77, 770) and 23 1(77, 770) vanish and, consequently, the asymptotic behavior of the 
beta functions is given by 23 2 (?7, 770) (assuming that it does not vanish). Hence, /3(?7, 770) ~ O (k~ 2 ), 
as we wanted to prove. Note that this result is valid for any value of fi and thus holds for every 
possible mass of the scalar field. In particular, it is so in the massless case, corresponding to 
ft = 3/2. 

Let us now study the beta coefficients /3k of the transformation that relates the Allen-Folacci 
vacua and the vacuum that defines the Fock representation with unitary evolution. Recalling the 
definition (|54l) of these beta coefficients, as well as the initial data that define both types of vacua 
[see expressions (l50l) . (152l) . and (|53l l. one arrives at 



/5k = -A k 



lojfc ■ 2 . 

■S(Xcos 770) - ^(cos 770) - i-S^cos r/o) - icoMcos 770) 



(BIO) 



n n 
As in previous calculations, we introduce the expressions of the functions Ry, Sy, and their time 
derivatives, in terms of the associated Legendre functions. Then, from the asymptotic expansions 
(14014 II) of those functions, and using Eq. (ISTI) . one finds that the beta coefficients admit an asymp- 
totic expansion of the form 



\ r(k -fi + 3/2) Y(k + fi + 3/2) v 

^~ \lnk + fi + 3/2) 2^ ix* +i) 2j^- (BU) 
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Here, 23 7 are functions of the asymptotic order O ~ O (k~^\ Employing again the Stirling 
formula (IB7I ), it follows that the global coefficient appearing in front of the sum of functions 23 7 is 
of order 0(1). Therefore, as in the previous case, the asymptotic order of the beta coefficients will 
coincide with that of the first nonvanishing contribution 23 ; . In this case, the function 23 reads 



23 n = cos 



(k + 1)770 + - + — 



sin 



n im 
(k+l)r l0 -- + Y 



+i sin 



(£+l)770 + - + y 



- I COS 



n n n fun 
(k+l)r]o-- + Y 



(Bl2) 



On the other hand, one can check that the function 33 1 is given by 
4/i 2 -! 



8(fc + 2) sin 7/o 
-i sin 



- cos 



+ sin 



n an 
(k + 2)7/0 + 4 + y 



(fc + 2)77 -- + ^- 



- I cos 



(* + 2)7/0 + ^ + y 



0. 



(Bl3) 



Therefore, as we wanted to show, the asymptotic behavior of the beta coefficients is at most of 
the order 23 2 ~ 0{k~ 2 ^. This guarantees their square summability (counting the degeneracy), and 
hence the unitary implementation of the transformation in the considered Fock representation. 
Again, this result is independent of the value of /j, and therefore is valid not only for the massless 
case {jjl = 3/2), but also for an arbitrary mass of the scalar field. 
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